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Abstract 


Computationally, selection of a proper numerical value for infinitesimal non Archime- 
dean epsilon in DEA models has some difficulties. Although there are several algo- 
rithms for selecting the proper non-Archimedean epsilon, it is important to introduce 
methods in order to calculate the efficiency of DMUs without using epsilon. One of 
these methods is a two-phase method, which obtains the efficiency of each DMU 
through solving two LPs, which the second LP is depended to the first. This paper 
proposes a method, which is able to compute the efficiency of DMUs by two LPs, 
which are not depended to each other and computationally can solve in a parallel 
computation. The major of this method is to find two references for each unit and 
combine them to obtain actual reference. 


Keywords: Data Envelopment Analysis (DEA); Decision Making Units 
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1 Introduction 


Since the first mathematical model of Data Envelopment Analysis (DEA) by 
Charnes et al. (1978), (known as CCR), and Banker et al. (1984) (known 
as BCC), there have been many theoretical and applied researches in DEA 
(Emrouznejad ea al., 2008). In 1979 the first version of DEA model has been 
updated by adding the non-Archimedean « as a lower bound for weights of 
inputs and outputs of the corresponding DMUs(Charnes et al. 1979). Differ- 
ent methods have been proposed for computing a suitable value for ¢. Ali 
and Seiford (1993) introduce a method to find an acceptable value for «. 
Mehrabian et al. (2000) modify this method and propose an LP to select a 
proper value for c«. Up to now, some researchers have published methods and 
discussions about the non-Archimedean ¢ such as Amin and Toloo (2004), 
MirHassani and Alirezaee (2004), Alirezaee and Khalili (2006). 

As the first and most important substitution method for epsilon-based 
DEA solving methods, Cooper et al. (1999) introduce the two-phase method 
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for evaluating the efficiency in DEA without using ¢. In this method , two 
LPs must be solved respectively for each DMU. 

In this paper, we introduce an alternative algorithm for evaluating the 
efficiency of DMUs without using ¢. In the proposed method, firstly, we find 
the references of DMUs and then inherited the references in a way that we can 
find out if the reference point of the unit is located on the weak frontier. The 
most advantage of this method is reducing the overall running time, because 
we can use parallel computation for our independent LPs in the algorithm. 


2 Classification of DMUs 


In DEA a set of DMUs are partitioned into two main classes: efficient and 
inefficient. The efficient units make the efficiency frontier. Figure 1 shows 
the classification of DMUs, based on the position of their reference on the 
efficiency frontier (Charnes et al., 1991). 


Fig. 1: Classifying DMUs. 


Tn this classification, EF and E’ are efficient and NE and NE are inefficient 
DMUs. In addition, F and NF are weakly efficient and weakly inefficient, 


respectively. 
Suppose that there are n Decision Making Units (DMUs) each consumes 
m inputs to produce s outputs. Let 2; = (%1j;,%2;,...,%mj) and y; = 


(Y1j; Ya7>-++sYszj) are input vector and output vector of DMUj (j =1,...,n), 
respectively. Hence, the CCR model corresponding to DMUp is as follow: 
CCR Model: 
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min z = 0 — e(7"4(s; ) + Vals, )) 
S.t. 
LipO — 5; _— ai Lig; = 0, Vi 


— st + Uri Ag = Yrp VT 
Aj,8, ,8¢ > 0,Vi,r, 9 


5134 9 Sp 


To introduce the new method, all of the values of variables for all DMUs are 
needed. So consider the following integrated model. 
CCRp Model: 


minz = 0, — e( L721 5;, + Pa Ae 


S.t. 
— = ig = Le LigAjp => 0, Vi 
i+ oF: “1YrjAjp =Urp, Vr 
ace roe Uy VOT 9 


In the above model, for each variable an index has been added. To simplify 
the notation, let u; = (y;,7;) be DMUj. Hence one can present members of 
productivity possibility set as u = (y, —a). Clearly, u; € PPS,(j =1,...,n). 
Let (05,A5. 93 ae | is an optimal solution of CC'R,, then efficiency of up 
is equal to 0. We also denote the efficiency of u by 6%. The dominate space 
and sefevedie of up, are denoted by DS, and u(p), respectively, and the set 


of reference indices of u, is denoted by E(p). 
yailj uj, E(p) = {J : Njp > OF 


It is clear that u(p) € DS, and u(p) = Yyex(p)Ajpuj- These concepts are 
illustrated in Figure 2. 
The shading pattern in this figure and other figure in the paper represents 
the dominant space of up, so every point in that region has inputs less than or 
equal and outputs more than or equal to u,. And dashed line in the frontier 
is weakly frontier. 

Therefore, we have: 


DS, = {u€ PPS|u > up}, u(p) = 


1. Up is efficient = 65 = 1, X72) 5;, e Pee = 0. 
2. Up is weak efficient = 6% = 1, 2 om ; se Sip >0. 
3. Up is inefficient = 05 <1 DP + a a a= = 
4, Up is weak inefficient = 05 < 1, Li 15;, a + Ds_jst a > 0. 


r= 


Notice that u(p) lies on the efficiency frontier, and so Cay. 
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Fig. 2: Efficient frontier, dominate space, reference point and set of reference 
index. 


3 New Method 


As explained in previous section, for classifying the DMUs, the values of 6 and 
slacks of reference point must be used. For each uy, two reference points are 
determined, the first reference belongs to DS, which minimize 6,. The second 
reference belongs to DS, which maximize sae + ana . Consider the 
following models: 


Model P1: 
min z = Ops.t. 
LipOn — bide B tigdby > 0, Vi 
a= ae = Urps Vr 
l 
Ajp 29, VG 
Model P2: 
min z = — (221%, + a ae) 


S.t. 
at is Liege = Lip, Wt 
t+ ai: aig AGe = =Yrp, Vi 


2 = 
2 ip? a0, Vj,0,7 
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The first reference for u, is determined by Model P1 as ui (P) = TEAL uy. 
Similarly, the second reference for u, is determined by the model p2 as 
u2(p) = BAG tie These references are illustrated in Figure 3: 


weakly efficient unit efficient unit 
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Fig. 3: up and its first and second references. 


u1(p) is not always the same as the reference point which is defined as 
(0*x;,y;) for (a;,y;). If up be an efficient or inefficient unit then its reference 
point and u;(p) are the same and if it is a weakly efficient or inefficient unit 
then probably its reference probably its reference point and u;(p) are not 
equal to each other. For example, in the Figure 4(a), u;(p) and reference 
point of u, are the same, but if the weakly efficient unit ui(p) is removed as 
illustrated in the Figure 4(b), the two definitions become different. In this 
case the reference point can be improved to the ui(p), therefore we must 
improve the reference point, too. Reference of each unit must be on the 
(strong) frontier and not on the weakly frontier. Weakly frontier is shown by 
dashed line in the figures in the paper. 


Definition 3.1. The revised reference for up is defined as bellow: 
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If u, be an efficient or an inefficient unit then the revised reference is u1(p), 
but if it is a weekly efficient or inefficient unit then the reference for up on 
the frontier is moved by i(p) using the weights belonging to E(P). Revised 
reference is simply the second reference of the first reference of the pth unit. 
We introduce a simple algorithm to identify the revised reference of up) and 
describe the idea of revised reference. It woks as follows: 


1. Find the first reference of u,. There are one or more units that create the 
first reference of u,. They are efficient or weakly efficient units. 

2. For each unit that participating the construction of the first reference of 
Up, find the second reference. For each of them, there are one or more units 
that are efficient (and not weakly efficient). 

3. The second reference of the first reference of up is the actual reference of 
it, so calculate values of the new variable that create the second reference 
of the first reference of u,. The revised reference is on the efficient frontier 
and is the actual reference of up. 

This algorithm computed the revised reference of uy, in other word we 
move of u, to ui(p) and then move to w2(wi(p)), as illustrated in Figure 
4. 
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Fig. 4: Reference point and ui(p) may be different. 


These concepts lead to the next theorems 1 proves that only efficient (not 
weakly efficient) units participating in construction of revised reference. And 
theorem 2 proves that comparing uy to the revised reference results the actual 
efficiency. 


Theorem 3.1. If uz is not an efficient unit then ey MAR =0, 
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Proof. Since rN; is the optimum weight of wu, in the model P2 for evaluating 
uj, if uz, is a weakly efficient, an —— or a weakly inefficient unit then 


for all j, Agj =0 and 7"_, ME Q2 = 


Theorem 3.2. If we select ti(p) as a reference for up in the model P1, then 
efficiency value of Up ts equal to 67. 


Proof. We knew that u;(p) is on the efficient frontier and 9% is the minimum 
value of 6,, so we must prove that t(p) > u1(p). We rewrite the u(p) as follow: 
a(p) = ey G, e324 Fux ] = a Appi 2(j), where u2(j) is the 
reference of wu; in the model P2. 
There are three possibilities for s;: 


1. If u; is an inefficient or a weakly inefficient unit, then Mo = 0. 
2. If u; is an inefficient unit, then ALD > 0 and ue(j) = uj. 
3. If u; is an weakly unit, then Me > 0 and ue(j) > uj. 


In all cases we have ti(p) = ye 1Ajpue(s) = as 1Ajpty = U(p). 
Since 7 is the minimum value of 6), therefore, if = select ti(p) as the 


reference of up in model P1, the efficiency value of u, is equal to 6). 


This shows that i(p) is on the efficient frontier, and )0/_, Anon > 0 if 
and only if u; is an efficient unit. 

Based on Theorems 1 and 2, t%(p) is a combination of (only) efficient units 
and its corresponding efficiency value is the same as in model Pl. After 
applying models P1 and P2 for all units, it is possible to compute the revised 
references of all units, which are on the efficient (and not on the weakly 
efficient) frontier. It sounds that this method is similar two phase method, 
but in the new method, the two models are independent. 

Based on theorems 1 and 2, t(p) is on the efficient frontier and as a ref- 
erence for u,, the efficiency value is not changed. a =F ii MAR where 
the related slacks are computed as follows: 


n n 

a= N¥e \ x af atx __ Ne oe 

Sip = Op Xip Dipti, Srp = pers Yrp- 
j=1 j=l 


Therefore, after applying Pl and P2 models for all units, we can compute 
the revised references, which are on the efficiency frontier, and the related 
slacks. According to classifying the DMUs, we have: 


1. If 0§ = 1 and 977", 85° + Deas Si = 0, then up is efficient, 

2. If or = = Land )yy", 8° + Do-1 Si > 0, then up is weakly efficient, 
3. If 6 < land yy") 8), + U1 Sty = 0, then up is inefficient, 

4. If 0 < land 7", 8° + Oo_, at* > 0, then wu, is weakly inefficient. 


1=1 “ip rp 
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The reference for up is i(p) = 054 Renee 

The basic role of this method is removing the non-Archimedean epsilon of 
models and using models that have no dependency to each others and could 
be solved separately while the models of two-phase method needs to be solved 
respectively, because second model uses of result of first model. 

The basis of traditional two-phase method is to find the optimum value 
of @ in the first stage and fix it to the second LP and solve it to find the 
maximum value for sum of slacks. 


4 Numerical Example 


In this section we solve a simple numerical example. We add the constraints 
jai jp = Land 0%_, AZ, = 1 to models P1 and P2 respectively and create 
the BCC versions of DEA models. 


Consider the following example: 


Table 1: Data for the numerical example 


Output Input 


y xz 
DMU1 1 1 
DMU2 2 1 
DMU3 _ 0.2 2 
DMU4 3 3 
DMU5 4 3 


These data are illustrated in Figure 5. 


Fig. 5: DMUs of the example. 
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Table 2 shows the optimum solutions of models Pl and P2 in the BCC 
format: 


Table 2: Results of the example using model P1 and P2 with variable returns 
to scale 


DMU | 6* dr 2" 
12345 [12345 
DMU1| 1 10000 /01000 
DMUvU2| 1 01000 /01000 
DMvu3/ 050/ 10000 101000 
DMU4)/0.6667/0 0.5 0 0 05/0000 1 
DMU5| 1 00001 0001 


For example for DMU4, we have \3, = 0.5, A2, = 0.5, and \2, = 1. The 
following table presents §~°, 8+” and A*: 


Table 3: 8~", +” and \* 


DMU 6* g— gt A 
T2345 
1 1 0.01001 000 
r= 000001 000 
3. 0.50 0.001800 1 000 
4 0.6667 0.0 0.0005 0 005 
5 1 000000 001 


Thus, the revised references of DMUs are as follow: 


Table 4: The revised references of DMUs 


DMU 1; u(3) EG) 
1 (,-l 1xw {2} 
2 (2,-1) 1x ug {2} 
3 (0.2, -2) 1X ug {2} 
4 (3,-3) 0.5 x ug +0.5 X us {2,5} 
5 (4-3) 1X us {5} 


It is concluded that we can compute the results of the efficiency evaluation 
of DMUs by applying a linear programming for each unit. 
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5 Conclusion 


There are two major methods for solving the basic DEA models: the epsilon 
based method, which selects a real number for epsilon, and the two phases 
method, which is used two LPs for each DMU. In this paper, we presented 
another method that determines two references for each DMU and then com- 
bines them and computes new lambdas and slack variables. Solving models 
without using non-Archimedean epsilon is an advantage of the new method 
and ability of computing the models in parallel can reduce the overall running 
time. 


References 


1. Ali, A.I. and Seiford, L. M., Computational Accuracy and Infinitesimals in Data En- 
velopment Analysis. INFOR, 37, (1993), 290-297. 

2. Alirezaee, M. R. and Khalili, M., Recognizing the efficiency, weak efficiency and inef- 
ficiency of DMUs with an epsilon independent linear program. Applied Mathematics 
and Computation, 183, (2006), 1323-1327. 

3. Amin, G. R. and Toloo, M., A polynomial-time algorithm for finding Epsilon in DEA 
models. Computers Operations Research, 31, (2004), 803-805. 

4. Banker, R. D., Charnes, A., and Cooper, W.W., Some models for estimating techni- 
cal and scale inefficiencies in data envelopment analysis, Management Science, 30, 
(1984), 1078-1092. 

5. Charnes, A., Cooper, W. and Rhodes, E., Measuring the Efficiency of Decision Making 
Units. European Journal of Operational Research, 2, (1978), 429-444. 

6. Charnes, A., Cooper, W. and Thrall, R. M., A Structure for Classifying and Char- 
acterizing Efficiency and Inefficiency in Data Envelopment Analysis. The Journal of 
Productivity Analysis, 2, (1991), 197-237. 

7. Charnes, A. and Rhodes, E., Short Communication Measuring the Efficiency of 
DMU’s. EJOR (1979), 339-339. 

8. Cooper, W. W., Seiford, L .M. and Tone, K., Data EnvelopmentAnalysis A Com- 
prehensive Text with Models., Applications, References and DEA-Solver Software, 
Springer Science, (1999). 

9. Emrouznejad, A., Parer, B. R. and Tavarese, G., Evaluation of research in efficiency 
and productivity: A survery and analysis of the first 30 years of scholarly literature 
in DEA. Socio-Economic Planning Sciences, 42, (2008), 151-157. 

10. Mehrabian, S., Jahanshahloo, G. R., Alirezaee, M. R. and Amin, G. R., An Assurance 
Interval for the Non-Archimedean Epsilon in DEA Models. Operations Research, 48, 
(2000), 344-347. 

11. MirHassani, S. A. and Alirezaee, M. R., An efficient approach for computing non- 
Archimedean (Epsilon) in DEA based on integrated models. Applied Mathematics 
and Computation, 166, (2004), 449-456 . 


